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Over the past five years, PHOBOS has collected data on proton-proton, deuteron- 
gold, Au+Au and Cu-|-Cu collisions, covering a wide range of collision energy, collision 
centrality and system size. Using these data, we have identified scaling features that are 
common for all types of high-energy collisions, as well as collective effects that are unique 
to the conditions created in collisions of relativistic nuclei. In this paper, we will focus on 
recent results obtained for collisions of Cu nuclei. Both in terms of universal features of 
particle production, and in the development of truly collective effects, the results for Cu 
nuclei confirm and extend our present understanding of nuclear collisions at the highest 
energies. In addition, we will describe recent unique results on multiplicity fiuctuations 
and particle production at very low transverse momenta. 

1. Overview 



PHOBOS is one of four experiments at the Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven National Lab, studying collisions of heavy nuclei at energies up to s/snn = 
200 GeV. Our goal is to study the properties of strongly interacting matter at extreme 
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temperature and density, where ab-initio numerical QCD calculations predict a phase 
transition to a system dominated by quark and gluon degrees of freedom. Our findings 
from the first four years of these studies have been summarized in [|T]. There we argued 
that in Au+Au collisions, a dense, highly interacting system is created, with energy 
densities far in excess of the critical values for the QCD phase transition. 

2. The PHOBOS experiment 

A detailed description of the PHOBOS detector can be found in The main com- 
ponents of the apparatus are a two-arm high-resolution magnetic spectrometer near mid- 
rapidity and a nearly A-ir charged particle multiplicity detector covering |?7| < 5.4. Both 
the spectrometer and multiplicity array use analog readout silicon pad sensors for particle 
detection, with a total of 135000 channels. The silicon detectors are complemented by 
scintillator based trigger and time-of-fiight counters, and forward calorimeters used for 
event characterization by detecting spectator nucleons. In five runs at RHIC, PHOBOS 
has collected more than 10^ events to tape. An overview of the available datasets can 
be found in Table ^ A key component of the experiment that is not usually discussed 
is the unique offline computing architecture that PHOBOS has established at the RHIC 
computing facility (RCF). We decided at an early stage to implement a distributed disk 
solution, equipping all nodes in our RCF Linux farm with cheap disk drives, for a to- 
tal storage capacity of close to 100 TByte. That allows us to keep our entire physics 
dataset available on disk for analysis at all times. In addition, we developed effective data 
management software and collaborated in the development of the PROOF extensions of 
ROOT [ISI, allowing us to use the main RCF cluster for interactive, parallel analysis of 
the large Cu-|-Cu and Au-|-Au data samples. The effective bandwidth achieved in these 
analyses using PROOF at RCF exceeds 5 GByte/sec. Details of the PHOBOS computing 
architecture can be found in [3i and a forthcoming publication. 
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Table 1 

Number of events stored to tape in millions for different systems from p+p to Au-|-Au 
and different collision energies ^/s^. 
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Figure 1. Ratio of the average number of col- 
lisions per participant, u, for Cu+Cu (light 
symbols) and Au+Au (dark symbols) as a 
function of Npart for y/s^ = 200 GeV from 
PHOBOS Glauber MC. 



Figure 2. Average eccentricity, estd, of the 
collision zone in Cu+Cu (light symbols) 
and Au+Au (dark symbols) function 
of Npart for = 200 GeV from PHO- 

BOS Glauber MC. 



3. Physics Results 

From the very beginning, the goal of PHOBOS was to obtain a broad survey of particle 
production in nuclear collisions for a large variety of system size, collision centralities and 
collision energies. As Table [T] shows, we have indeed collected the datasets to fulfill this 
original program. Using these data, we attempt to extract organizing principles or scaling 
features directly from a systematic study of the available observables. As we will show 
in the following, several striking scaling rules have indeed emerged from the data. Some 
of those appear to connect the results obtained in nucleus-nucleus collisions with those in 
more elementary collisions, like p+p, p+A or e~^e~. Other results show the existence of 
unique collective effects not present in elementary systems. Clearly, the true significance 
of these features will be established when they can be shown to naturally emerge from 
the eventual theoretical understanding of the physics of these collisions, or when future 
experiments at much higher energy at the LHC can be shown to be constrained by the 
same underlying features as data at present energies. 

In the following subsections, we will use the PHOBOS data to address three distinct 
questions about scaling of particle multiplicities and distributions in nucleus-nucleus colli- 
sions. The first question concerns the scaling of overall particle production in Cu+Cu and 
Au+Au collisions. In the following section, we will investigate the connection between 
the centrality and energy dependence of particle production and finally we will study the 
connection between of the final-state azimuthal anisotropy and the initial state geometry. 

As we will show in the following, conclusive answers emerge from the data for all these 
questions. It is important however to consider what differences between Au+Au and 
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Figure 3. Charged hadron pseudorapidity distributions for different centralities for 
Au+Au collisions at y^s^ = 19.6, 62.4, 130 and 200 GeV (top row from left to right), 
Cu+Cu at 62.4 and 200 GeV (middle row, preliminary) and d+Au at 200 GeV (bottom 
row). 



Cu+Cu could be expected based on the underlying differences in collision geometry. For 
this, the result of our Glauber calculations plotted in Figs. [T] and El are important. Shown 
in Fig. ^ is the average number of binary NN collisions u every nucleon undergoes within 
a nucleus-nucleus collision, as a function of the number of participating nucleons Npart- 
As this plot shows, a study of the centrality dependence of particle production allows a 
large variation of u. However, for a given value of Npart, ^ in Cu+Cu and Au+Au is the 
same within the possible experimental resolution and therefore the Cu+Cu vs. Au+Au 
comparison does not provide additional discriminating power for distinguishing Ncoii vs. 
Npart scaling. To better discriminate effects scaling with the number of binary collisions 
Ncoii from those scaling with Npart, a comparison of Au+Au results with those from 
much smaller nuclei like e.g. silicon or carbon would be necessary. However, Fig. |21 shows 
that there is a very important difference between the collision geometry in Cu+Cu and 
Au+Au, even at the same Npart- In this figure, the average eccentricity, estd, of the collision 
zone in the transverse plane is plotted against Npart- Naturally, for the same Npart the 
eccentricity in peripheral Au+Au collisions is much larger than in the corresponding more 
central Cu+Cu collisions. The comparison of Cu+Cu and Au+Au will therefore allow a 
direct investigation of the connection between initial density and the efficiency with which 
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Figure 4. Pseudorapidity distribution for 
charged hadrons in Cu+Cu collisions (black 
symbols) and Au+Au collisions (open sym- 
bols) at = 62.4 GeV. The Cu+Cu 
and Au+Au centralities were selected to 
yield similar Npart- The grey band indicates 
the systematic uncertainty for Cu+Cu (90% 
C.L.). Errors for Au+Au are not shown. 



Figure 5. Pseudorapidity distribution 
for charged hadrons in Cu+Cu collisions 
(black symbols) and Au+Au collisions 
(open symbols) at y/s^ = 200 GeV. 
The Cu+Cu and Au+Au centralities were 
selected to yield similar Npart- The 
grey band indicates the systematic uncer- 
tainty for Cu+Cu (90% C.L.). Errors for 
Au+Au are not shown. 



the collision system translates the initial state geometrical asymmetry, characterized by 
Cstd, into a final state asymmetry in the momemtum distributions. The results of this 
comparison, which require a detailed understanding of the definition of the appropriate 
eccentricity measure, will be discussed in detail in section 13.31 

3.1. System-size dependence of particle production 

One of the main questions for the Cu+Cu run at RHIC in early 2005 was that of 
the system-size dependence of particle production in nuclear collision at RHIC energies, 
in terms of overall multiplicity, dN/drj distributions and dN/dpx distributions. Before 
discussing in detail some of the scaling features observed in the data. Fig. El gives an 
impression of the breadth and quality of the data collected by PHOBOS. The figure 
shows charged hadron pseudorapidity distributions as a function of centrality for Au+Au 
collisions at = 19.6, 62.4, 130 and 200 GeV [ElEl, Cu+Cu at 62.4 and 200 GeV 

and d+Au at 200 GeV [|7|. 

For both Au+Au and Cu+Cu collisions, the centrality selection was based on the 
distribution of energy deposited in the region of the multiplicity counter covering |?7| < 3.0 
and in the trigger counters covering 3.2 < \ri\ < 4.5. The energy distributions for a 
minimum bias dataset are divided into bins of fractional cross-section and the average 
number of participants for each bin is determined from a Glauber calculation using the 
HIJING event generator [|H] and a GEANT-based simulation of the PHOBOS detector 
response. Details of this method can be found in [ [T] . 
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Figure 6. Nuclear modification factor Raa as a function of Npart in bins of px for collisions 
aty'Sjvjy = 200 GeV near mid-rapdity. Data for Cu+Cu collisions (filled symbols) are 
shown in comparison with data for Au+Au collisions from PHOBOS (open circles) and 
PHENIX (open diamonds). Systematic errors for Cu+Cu data (90% C.L.) are shown as 
brackets. 



For the comparison of dN/drj distributions in Cu+Cu and Au+Au, we selected cen- 
trality bins in Cu+Cu and Au+Au such that the average number of participants approx- 
imately matched. The result of this comparison is shown in Figs. 0] and El for collisions 
at 62.4 GeV and 200 GeV, respectively. No further scaling factors are applied. This 
comparison, and further studies for different centrality bins, provide a strikingly simple 
answer to the question above: If Cu+Cu and Au+Au collisions at the same collision 
energy are selected to have the same (Npart), the resulting charged particle dN/drj distri- 
butions are nearly identical, both in the mid-rapidity particle density and the width of 
the distribution. This is true for both 62.4 GeV and 200 GeV data. 

Following the similarity in both the mid-rapidity density and the shape of the pseudora- 
pidity distribution, it is natural to compare particle production as a function of transverse 
momentum in Cu+Cu and Au+Au, for different centrality classes. This comparison can 
be seen in FigEl where we plot the nuclear modification factor Raa as a function of Npart 
for Cu+Cu and Au+Au collisions at 200 GeV, in bins of from 0.25 to 6.25 GeV/c. 
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Raa is defined as 



^AA{PT = TTT Tl 1— • [i-) 



For Au+Au, data from PHOBOS [ E] and PHENIX [ CHI are sfiown. It is important to 
note again tliat tlie same Npart value in Cu+Cu and Au+Au corresponds to virtually 
the same value of Ncou for both systems. Therefore, the fact that Raa coincides for 
Cu+Cu and Au+Au at the same Npart, as shown in Fig. IHl implies that for a given system 
size, measured by either Npart or Ncou, the absolute yield per participant is the same in 
both systems. This statement holds over the entire range in centrality and pt covered 
in this analysis and is also confirmed with similar precision in collisions at 62.4 GeV. In 
summary, overall particle production per participant, as well as particle distribution in 
pseudorapidity and transverse momentum, appear identical in Cu+Cu and Au+Au at 
a given collision energy, if collisions with the same average number of participants are 
selected. 

3.2. Factorization of energy and centrality dependence 

While the comparison of particle production in Cu+Cu and Au+Au at the same col- 
lision energy yields the simplest possible connection between the two systems, a simul- 
taneous study of centrality and energy dependence of the same observables for either 
of the systems yields a subtle and surprising scaling relationship, that again holds for 
mid-rapidity particle production, the shape of the dN/drj distributions and the transverse 
momentum distributions. As was first shown for the mid-rapidity density [E], the ob- 
served increase in particle production per participant with increasing Npart (see FiglTj) is 
independent of collision energy over the full energy range of RHIC from 20 to 200 GeV. 

This is illustrated in FiglSJ where we show the ratio of mid-rapidity densities as a func- 
tion of Npart relative to 200 GeV data. All ratios, from 200/130 GeV to 200/19.6 GeV, 
are fiat within the experimental uncertainty. This factorization of energy and centrality 
dependence is strongly violated in models that describe overall particle production as a 
superposition of independent "soft" contributions scaling like Npart and "hard" contri- 
butions scaling like Ncou, where the increase in the hard contributions is given by the 
increasing mini-jet production with energy. As a quantitative example. Fig. |H1 shows 
the prediction from HIJING. It is interesting to note that the observed factorization is 
described in approaches based on the ideas of parton saturation [ IT^ ITHj. 

3.2.1. Energy and centrality factorization of pt distributions 

In an earlier publication [|^, we showed that the factorization described above not 
only holds for the pT-integrated yields at mid-rapidity, but also differentially over a large 
range of transverse momenta. This is illustrated in the top row of Fig. El for Au+Au 
collisions at 62.4 and 200 GeV. There, the quantity Rpq'''\pt), defined as 

, _ iK-f) d'NAA/dpTdy 

^"^^^ ~ {Npart) dmr/^ydprdv ^ ^ 

is shown as a function of px for the six centrality bins. Remarkably, the centrality evolution 
of the spectral shape from central to peripheral collisions is found to agree within the 
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Figure 7. Mid-rapidity density 

dN / drj/ {Npart/'^) as a function of Npart- 
Au+Au data for = 19.6, 62.4, 130 

and 200 GeV are shown as filled symbols, 
Cu+Cu data for ^/s^ = 62.4, and 200 GeV 
are shown as open symbols. Ellipses indi- 
cate the uncertainty in Npart and dN/drj 
determination. Lines indicate the additional 
uncertainty due to the uncertainty in the 
trigger cross-section for Cu+Cu. 



Figure 8. Ratio of mid-rapidity den- 
sities dN/dr]/{Npart/'2) function of 
Npart- Au+Au data for energies 200/19.6, 
200/62.4, 200/130 (filled symbols) and 
Cu+Cu data for 200/62.4 (open symbols) 
are shown, in comparison to HIJING cal- 
culations and predictions from Kharzeev 
et al. for the 200/19.6 ratio [ 12 . Lines 
indicate the additional uncertainty due to 
the uncertainty in the trigger cross-section 
for Cu+Cu. Ratios for p+p collisions are 
shown at Npart = 2. 



experimental uncertainty between 62.4 and 200 GeV, even though Raa itself shows a large 
variation by more than a factor of two at high pt between the two energies. The bottom 
row of Fig. El shows that the same energy-independence of the shape evolution is seen for 
the much smaller Cu+Cu system, when again plotting Rpcipr) in bins of fractional cross- 
section. This apparent dominance of collision geometry, relative to expected dynamical 
effects, is an important feature of the data that remains to be understood. 

3.2.2. Extended longitudinal scaling 

As shown in earlier publications, the energy and centrality dependence of particle pro- 
duction factorize not only for mid-rapidity multiplicity densities. One of the most striking 
scaling relationships seen in nuclear collisions at RHIC is the energy-independence of par- 
ticle yields at moderate to high rapidities, when viewed in the restframe of one of the 
colliding nuclei. This phenomenon has also been observed in p+p and p+A collisions 
over a large range of energies and is commonly referred to as "limiting fragmentation" . 




Figure 9. Ratios of yields per participant vs pt (near mid-rapidity), relative to central 
collisions for bins of fractional cross- sect ion. Central events are shown in the right-most 
panel. Solid symbols show 200 GeV data, open symbols show 62.4 GeV data. Top row 
shows Au-|-Au data, lower row shows Cu+Cu data. Systematic uncertainties (90% C.L.) 
are shown by brackets. 



PHOBOS has extended this observation not only to dN/drj distributions in nuclear colli- 
sions at different colhsion centralities [IH], but also to the pseudorapidity dependence of 
azimuthal anisotropies [^mHl- In Fig. a compilation of results for dN/drj, the elliptic 
flow coefficient V2 and the directed flow coefficient vi is shown for Au+Au collisions at 
four different collision energies from 20 to 200 GeV. 

As can be seen in this figure, limiting fragmentation holds not only for the multiplicity 
distribution of charged hadrons, but also for observables like V2- This is remarkable as 
V2 is believed to be generated dynamically in the evolution of the collision, from pre- 
sumably very different initial conditions at the different collision energies. Furthermore, 
Fig. ITUl shows that the scaling behaviour is not confined to a small "fragmentation region" 
near beam-rapidity, but extends over a large part of longitudinal phase space. For this 
reason, we refer to this set of observations as "extended longitudinal scaling". Finally, 
Fig. im shows that the same scaling behaviour also holds for dN/drj and V2 in Cu+Cu. 
Measurements of vi in Cu+Cu are forthcoming. 

In order to examine the connection between energy- and centrality-scaling for longitu- 
dinal distributions, it is useful to again use the ratio of yields for peripheral relative to 
central events, RpQ''\r]): 



,iv,„..,^, _ (A^°a-.f°) dNAA/dy 

{Npart) dNl-f^/dT] 



R%-iv) = vr / r:;-:. : (3) 



In Fig. ^1 we show that also the longitudinal distributions, when viewed in the restframe 
of one of the colliding nuclei, exhibit the factorization seen for mid-rapidity yields. 

In summary, the combined study of energy- and system size dependence of hadronic 
observables in Au+Au collisions has revealed a remarkable factorization, leading to an 
energy-independence of the centrality evolution for the mid-rapidity particle density, the 
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Figure 10. Dependence of dN/drj' (top), el- 
liptic flow t>2 (middle) and directed flow Vi 
(bottom) on r/' = \ri\ — yheam for Au-fAu 
collisions at ^fs^ = 19.6, 62.4, 130 and 
200 GeV. 



Figure 11. Dependence of dN/drj' 
(top) and elliptic flow V2 (bottom) on 
^' = I'?! ~ Vbeam for Cu+Cu coUisious at 
yJ^NN — 62.4 and 200 GeV (preliminary). 



charged hadron px distributions near mid-rapidity and the extended longitudinal scaling 
of multiplicity distributions. This factorization is particularly surprising in the context 
of a hard/soft picture of particle production, where the increasing mini-jet cross-sections 
with energy and the increasing ratio of collisions per participant, z/, couple energy and 
centrality dependence. It appears that the observed dominance of collision geometry, 
which is naturally energy-independent for a given centrality selection, is better captured 
in a description of the collision process based on ideas of parton saturation. 

3.3. System size dependence of elliptic flow 

The final question to be answered by the comparison of Au+Au and Cu+Cu data as a 
function of energy and centrality concerns elliptic flow, and the connection between the 
initial state conditions and the observed final state anisotropy. It has been argued in the 
past [HZ], that the observed flow coefficient V2 for Au+Au closely follows the initial state 
eccentricity as a function of centrality. This can be seen by comparing Npart dependence 
for the Au+Au initial state eccentricity from a Glauber calculation shown in Fig.|21and the 
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Figure 12. Ratio Rpcirf) of dN/drj' distri- 
butions for peripheral to central events for 
Au+Au collisions at y/s^ = 19.6, 130 and 
200 GeV as a function of 77' = \ri\ — ybeam- 
Only systematic uncertainties (90% CL.) are 
shown. 



Figure 13. Average elliptic flow coefficient 
V2 measured near mid-rapidity, func- 
tion of Npart for Cu+Cu collisions (open 
symbols) and Au+Au collisions (filled 
symbols) at ^/s^ = 200 GeV. Only sta- 
tistical errors are shown. 



Au+Au f 2 coefficient at mid-rapidity shown in Fig. ^| However, doubt on this crucial 
connection is cast by comparing the Cu+Cu calculation and data in the same figures. 
Whereas the average eccentricity for Cu+Cu tends to zero for the most central collisions, 
the corresponding V2 values for Cu+Cu only drop to a value of V2 ~ 0.03 even for the 
most central collisions. This would lead to the paradoxical conclusion that for the same 
Npart and therefore the same initial area density of produced particles, the Cu+Cu system 
is much more effective in translating an initial eccentricity into a final state anisotropy 
than the Au+Au system. Alternatively, large non-flow effects mimicking a dynamically 
generated anisotropy could be postulated for the Cu+Cu system. 

However, a possible explanation unifying the observations in Au+Au and Cu+Cu can 
be found by examining the underlying deflnition of the initial state eccentricity. This 
eccentricity, called egtd in the following, is commonly deflned as the average eccentricity of 
the distribution of participating nucleons, relative to the known reaction plane, obtained 
for a certain centrality class in a Glauber calculation. This deflnition suffers from two 
potential problems: It averages out the fluctuations from event-to-event in the actual 
participant distributions. Finite number fluctuations will lead to an eccentric nucleon 
distribution even for collisions with impact parameter 6 = 0. With the standard deflnition, 
these fluctuations will be averaged to zero for central events. Furthermore, the minor axis 
for the actual event-by-event participant distribution will in general not coincide with the 
impact parameter vector. The eccentricity calculated relative to the reaction plane will 
therefore underestimate the true eccentricity of the nucleon distribution. To study these 
deflciencies, we have deflned an alternative measure of the eccentricty in each centrality 
bin, where we calculate the eccentricity for each Glauber event relative to the principal 
axes of the actual participant distribution (see [d])- By construction, this participant 
eccentricity , epart, will always be positive and will therefore average to a flnite value 
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Figure 14. Participant eccentricity epart of 
the collision zone in Cu+Cu (light symbols) 
and Au+Au (dark symbols) function of 
Npart for = 200 GeV from PHOBOS 

Glauber MC. 



Figure 15. Ratio of V2 coefficient 
to average participant eccentricity as 
a function of mid-rapidity density 
l/{S){dN/dy) for Au+Au collisions at 
^/s^ = 130 and 200 GeV and Cu+Cu 
collisions at y/s^ = 62.4 and 200 GeV, 
in comparison to lower energy data from 
NA49 and E877. 



even for the most central events. In addition, the smaller number of colliding nucleons, 
makes the difference between egtd and epart particularly important for the Cu+Cu system 
relative to Au+Au. The result of a Glauber calculation for epart for Cu+Cu and Au+Au 
as a function of Npart can be seen in Fig. HH As expected, epart remains finite even for 
the most central Cu+Cu collisions. 

Using epart, we can now attempt to identify a common scaling behaviour of Cu+Cu and 
Au+Au collisions over a large range of collision energies and centralities. This is shown in 
Fig. plotting the ratio of (f 2) / (epart) versus the mid-rapidity area density of produced 
particles [021 1201 ■ The data appear to exhibit a common scaling behavior over a large 
range in collision energy, suggesting that the efficiency for translating the initial state 
eccentricity estimated using epart into a final state anisotropy f 2 appears to only depend 
on the initial area density achieved in the collision. Clearly, it is a fascinating question for 
future experiments whether this curve saturates at higher densities or continues to rise. 

4. Future PHOBOS physics program 

The future physics program of PHOBOS will rely on the analysis of the existing 
datasets. The program is characterized by two major directions. We will use the unique 
acceptance of our mid-rapidity spectrometer to study particle production at the lowest 
transverse momenta accessible at RHIC, as a function of collision energy, centrality and 
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colliding species. This includes studies of 0-meson production at very low px, which re- 
quires an analysis of the full high statistics datasets for Cu+Cu and Au+Au. The second 
class of studies uses the high statistics datasets in combination with the large acceptance 
of the multiplicity array to study fluctuations and correlations in particle production. 
Preliminary results from these studies are presented in separate contributions to these 
proceedings [ |^ 122 123 and will be summarized here briefly. 

4.1. Particle production at low px 

The unique setup of the PHOBOS silicon spectrometer and its proximity to the interac- 
tion vertex allows us to detect charged particles that range out and stop in the front layers 
of the silicon detector. For the present analysis, we have developed an algorithm detecting 
particles that stop in the fifth silicon layer. Based on the total deposited energy and the 
specific energy loss in the silicon detectors, we can reconstruct the momentum and mass 
of charged particles. This allows us to measure pions, kaons and protons (summed with 
their anti-particles) for transverse momenta as low as 30, 90 and 140 MeV/c, respectively. 
For details of the reconstruction algorithm, see [ I24j . 

At this conference, new results for the centrality dependence of particle production at 
low pt in Au+Au collisions at 62.4 and 200 GeV were presented, as well as results for 
minimum bias d+Au collisions at 200 GeV. A detailed discussion of the results can be 
found in [|^. As an example. Fig. [TBI shows the yield for pions, kaons and protons at 
low Pt in combination with results at higher pt for minimum bias d+Au collisions at 
200 GeV. Both results are fully corrected, including a correction for feed-down from weak 
decays. Also shown are the results of a blastwave fit to the higher pj- data, which is 
found to underestimate the low px points for d+Au, whereas for similar fits in Au+Au 
collisions consistency between the low and high pt datasets is found. It is interesting to 
note however that the effective expansion velocity extracted from the fit to the d+Au 
high Pt data is 0.35 c, indicating that the physical interpretation of the extracted value 
has to be treated carefully. 

4.2. Fluctuations and Correlations 

The high statistics Au+Au dataset, in combination with the large PHOBOS acceptance, 
allows us to perform a search for events that in a statistically quantifiable way differ from 
average Au+Au events. Possible mechanisms for the occurrence of such events might be 
the formation of droplets due to supercooling [123 or the formation of disoriented chiral 
condensates [123- While the likelihood of such scenarios is unclear, an unbiased search for 
"unusual" events clearly is an important part of the RHIC physics program. The details 
of our strategy for identifying rare events is described in [123- Using the most central two 
million events of our Au+Au dataset, we determine the shape of the average uncorrected 
dN/drj distribution for events in fine bins of vertex position. Similarly, the variance 
around the average shape in bins of t] is extracted from the data. Using the average shape 
and variance obtained from the data, we then calculate the of each individual event 
relative to the ensemble average. The resulting distribution is shown in Fig. El The 
distribution consists of a central core close to a x^/DoF ^ 1, with a tail out to large 
X^. In total, about 0.01% of all events are found in a region where for purely statistical 
fluctuations no entries would be expected. Further studies of these unusual events have 
shown that their rate is linearly related to the instantaneous luminosity at which the 
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Figure 16. Transverse momentum spectra 
for pions, kaons and protons in 200 GeV 
min-bias d+Au collisions. Lines indicate the 
results of a blast-wave fit for data above 
Pt = 0.25 GeV/c. 



Figure 17. Reduced distribution 
for event-by-event comparison of central 
Au+Au collisions at 200 GeV to ensemble 
average. 



collisions were recorded. At present, we therefore have no evidence for the existence of 
unusual physical fluctuations. Further studies are underway to set quantitative limits on 
various physical scenarios. 

4.3. Forward-backward multiplicity correlations 

Further examination of Fig. fTTI shows that the width of the distribution around unity 
is significantly larger than expected for a purely Poissonian production and detection of 
charged particles. Thus particles appear to be produced in a correlated fashion, rather 
than one-by-one. This can be studied quantitatively using correlations of multiplicities in 
non-overlapping bins of pseudorapidity. We have performed such studies for symmetric 
pairs of pseudorapidity bins centered at values between 0.25 < r] < 2.75, varying the bin 
width from 0.5 < At] < 2.0. To quantify the relative fluctuations between the multiplicties 
Np in the forward [rj > 0) and Nb in the backward [r] < 0) bins, we define the event- 
wise observable C = {Np — Nb)/ ^/Np + Nb- This variable has the useful property that 
its variance is one for independent particle emission, even when averaged over events 
from centrality bins of finite widths. After correction for detector and acceptance effects, 
which are described in [ EH 122] , "^^^ t)e used to study short range correlations in 
particle production. If particles are produced as clusters which decay with a rapidity 
width smaller than the typical bin width chosen in the analysis, then a"^ in the absence 
of other correlations will directly correspond to the cluster size /c, i.e. the multiplicity of 
decay products from each cluster. In Figs.^land^l we show the dependence of on the 
position and width of the bins used in our analysis for central Au+Au events. The main 
result is that cr^ is much larger than unity, in particular for larger A?], indicating that 
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Figure 18. Dependence of the fluctuation 
measure o"^ on the position rj for Arj = 
0.5 of the forward and backward multiphc- 
ity bins for central Au+Au collisions at 200 
GeV (preliminary). Systematic uncertainties 
(90% C.L.) are shown as brackets. 



Figure 19. Dependence of the fluctuation 
measure on the width At] for t] = 2.0 
of the forward and backward multiplicity 
bins for central Au+Au collisions at 200 
GeV (preliminary). Systematic uncertain- 
ties (90% C.L.) are shown as brackets. 



particles in Au+Au collisions are not produced independently, but rather in clusters. The 
results are reminiscent of those obtained in similar analyses for p + p collisions [EHj- The 
similarity of the results in A+A and p+p collisions, as well as the weak energy dependence 
seen in p+p, could indicate that the cluster formation is a phenomenon related to common 
features at hadronization. This will be further tested by future studies of multiplicity 
correlations in Au+Au and Cu+Cu collisions as a function of collision energy. For further 
discussion and a comparison of the present results with event generators, see [1^. 



5. Summary 

In this paper, we have presented new results from the recent Cu+Cu run at RHIC. The 
data show that particle production per participant nucleon is very similar in Cu+Cu and 
Au+Au collisions, if one selects collisions with equivalent number of participants. This is 
true not only for the mid-rapidity particle density, but also for px spectra and the shape of 
the pseudorapidity distributions. Like Au+Au collisions at RHIC, the Cu+Cu data also 
exhibit several intriguing scaling relationships, including the extended longitudinal scal- 
ing of pseudorapidity distributions and elliptic flow, and the factorization of energy and 
centrality dependence. The flrst data on elliptic flow in Cu+Cu reveal a shape of ^2(^7) 
that is similar to observations in Au+Au, with a surprisingly large magnitude of V2 near 
mid-rapidity relative to the expected average initial state anisotropy. We have argued 
that the result can be quantitavily understood, when taking fluctuations in the initial 
state transverse geometry into account. This leads to a universal scaling of V2 relative to 
the properly deflned "participant eccentricity" in Cu+Cu and Au+Au data over a large 
range in collision energy. We also presented examples of the ongoing systematic study 
of particle production at very low transverse momenta. For Au+Au, the low pt data 
agree well with extrapolations from higher pt assuming a radially expanding source. We 
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expect further insight into the underlying mechanisms of particle production from mea- 
surements of event-by-event fluctuations. Results were shown that revealed significant, 
cluster-like, correlations in final state multiplicity production, but no indication for the 
presence of unusual, large scale fluctuations. Overall, our data point to the importance 
of understanding the dynamics of the very early stage of the collision, where the scaling 
features observed in the final state hadron production and anisoptropies are established. 
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